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■ The self-energy and vertex QED radiative corrections to the parity non- 

conservation (PNC) amplitude in atoms are obtained using the perturbation 
theory in powers of otZ. The calculated linear in aZ term gives —0.6% for 
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the PNC amplitude in Cs. The estimated nonlinear terms make corrections 
larger —0.9(2)%. This result brings the experimental data for the 6s — 7s 
transition in 133 Cs in agreement with the standard model. 



> 1 32.80.Ys, 11.30.Er, 31.30.Jv 

00 

o 
o 

a: 
hi. 



X 



Typeset using REVTeX 



* kuchiev@newt . phys .unsw.edu.au 



1 



This work is inspired by a recently discovered pronounced deviation of experimental 
data on parity nonconservation (PNC) in atoms from predictions of the standard model. 
The deviation was found in the 6s — 7s PNC amplitude in 133 Cs, that was investigated by 
Bouchiat and Bouchiat Gilbert and Wieman ||, and, most recently with the highest 
precision 0.3%, by Wood et al ||. The experimental progress was matched by atomic PNC 
calculations that are crucial for the analysis of the experimental data. Accurate calculations 
of Refs. have been revisited recently by Kozlov et al || and Dzuba et al Bennett 
and Wieman || compared the theoretical data with available experimental data on 
dipole amplitudes, polarizabilities and hyperfme constants for Cs, and suggested that the 
theoretical error for the PNC amplitude should be reduced from 1% to 0.4%. This suggestion 
is supported by Ref. §| that estimates the theoretical error as 0.5%. There have been 
discovered recently several, previously neglected phenomena that contribute at this level of 
accuracy. The Breit corrections give —0.6% as was found by Derevianko |10| and confirmed 
in 



[T||. Sushkov |T2|] proposed that the radiative corrections may be important. This 
suggestion is supported by the calculations of Johnson et al |L3| who demonstrated that the 
QED vacuum polarization gives 0.4%, the value confirmed in flll||?|,|l~5|| . Ref. indicates that 
there is a 2.3a deviation of the weak charge Qw extracted from the atomic PNC amplitude 
from predictions of the standard model |2"Bfl . More recent analyses in which the 



Breit corrections (—0.6%) and the QED vacuum polarization (0.4%) were included, give 
similar deviations 2.2a and 2.0a respectively. 

We evaluate the self-energy and vertex radiative corrections, called the e-line corrections 
below. They were considered previously by Marciano and Sirlin and Lynn and San- 



dars using the plane wave approximation that resulted in a small value ~ 0.1%. An 



attempt to take into account the strong Coulomb field in this problem was made in Ref. [14 
that concluded that the e-line corrections are negligible because they are not enhanced by 
large logarithmic factors (which play an important role for the vacuum polarization, as was 
demonstrated by Milstein and Sushkov in the cited paper and confirmed in Ref. |T3[). We 
show that the strong Coulomb field makes the e-line corrections large due to unexpected 
numerical interplay. Our results are in qualitative agreement with conclusions of recent 
Ref. [|18j that demonstrated for the first time that the e-line corrections are important. We 
find quantitative agreement with this paper as well, the correction —0.9% obtained for the 
6s — 7s PNC amplitude in 133 Cs in the present paper matches the result —0.7(2)% of 



This correction brings the experimental data of Wood et al in agreement with the stan- 
dard model. The present paper relies on the perturbation theory in powers of aZ. It should 
be noted that |1| uses very different approach, based on an identity that expresses the 



radiative corrections to the PNC amplitude via corrections to the energy shifts due to the 



finite nuclear size that have been calculated previously in Refs. fl20|-f2~3|j . 

Diagrams (1),(2) and (3) in Fig.l describe the e-line radiative corrections to the PNC 
amplitude in the first order of the perturbation theory in aZ . The weak PNC interaction 
between the electron and the nucleus originates from the Z-boson exchange. It is described 
by the Hamiltonian (relativistic units h — c — m — 1 are used, if not stated otherwise) 

tfp NC = (2V2)- 1 GfQw p(r) 75 , (1) 

where Gf is the Fermi constant, Qw is the nuclear weak charge, and the nuclear density for 
our purposes can be approximated by the delta function, p(r) = Z5(r). The perturbation 
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induced by the Hamiltonian ([[]) is described by a factor GpQ^Zj (2\/2)7o75 that appears for 
all the diagrams considered. Since we are interested merely in a relative contribution of the 
corrections, the constant GpQ^Z '/(2v2) will be dropped in all formulas below to simplify 
notation. We use the Feynman gauge for the photon propagator D^ u (q) = Ang^/ (q 2 — A 2 ), 
and introduce the density matrix p = (1/4) (p + 1)(1 + 0,75) that distinguishes the electron 
state with the given momentum p^ and the polarization four- vector a^. The latter can be 
expressed via the polarization three-vector £, a = p • £, a = £ + p(p • C)/(l + Po)- It suffices 
to discuss the forward scattering in which the electron momentum p^ and its polarization 
remain intact, and restrict our consideration to the limit p — > because all important events 
that contribute to the process happen at separations r ~ 1/m that are much less than the 
atomic radius. The only pseudoscalar available in this kinematics is ao = p • C that should 
appear as a common factor for all diagrams for the PNC amplitude. We need therefore to 
extract the lowest, linear in p term of the p — > expansion. The analytical expression for 
each Feynman diagram can be written in the following form 

. f d 4 q d 3 k 4tt -AnZa Tr(A/» 
~ l0t J (27r) 4 (27r) 3 q 2 - A 2 ~k 2 V ' ^' 

where q is the photon momentum, k is the momentum running over the Coulomb and the 
PNC weak interactions, Tr is the trace over the spinor indexes, and the nominator M and 
denominator T> are specific for each diagram. To clarify notation let us present explicitly 
the nominators Mi, Mi for one Feynman diagram of the type (1) and one diagram of the 
type (2) in Fig.l 

Mi = Y(P + 1)70(6 + 1)7„(£ + l)7o75 , (3) 
Mi = Y(P + l) 7 o(Q + l)7o7 5 (P + 1)7/, , 

where P = p + q, Q = p + k + q, R = p + k, k — (0, k). The vertexes 70 and 7075 in Eq.@ 
originate from the Coulomb and the PNC weak interaction respectively. Using Eqs.( |Al| ) we 
find from Eq.(Q) (and similar expressions for the Feynman diagrams that are topologically 
identical to either diagram (1), or (2) in Fig.l) 

Tr(JVip) = " 8 [ 2P o(«o + (pR)a - (aR)p ) - (QPRa) ] , (4) 
J2 Tr(A/" 2 p) = -4 [ (Qa) - (PQRa) } . (5) 

Here (xy) is a scalar product of vectors, (xy) = x^, and the symbol (xyzu) stands for the 
following invariant of four four- vectors (xyzu) = (xy)(zu) — (xz)(yu) + (xu)(yz). The tilde 
sign in Eq.(^) marks the inversion of a four- vector x M = (xq, — x). Summation in (^) and 
(j|) includes all Feynman diagrams topologically identical to diagram (1) and (2) in Fig.l 
respectively. The conventional parameterization permits the denominator for diagram (1) 
in Fig.l to be presented in the form 

D 1 = ei(P a -i)+6(g 2 -i)+e8(9 a -A a ) . 
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Using similar presentation for diagram (2) we calculate the integral over the photon mo- 
mentum q in Eq.(|2|) for these two diagrams with the help of Eqs.(A2). Expanding the found 
results up to the first power in p we integrate them over the angular variables 1\ of the 
momentum k. These straightforward calculations lead to the following representation for 
the e-line radiative correction <5™^ e to the PNC amplitude 



PNC 
e— line 
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^PNC 



C PNC a 2 Z 



37T 2 Jo 



fi(k 2 ) + f 2 (k 2 ) + h(k 2 )) dk . 



(7) 
(8) 



Note that (5™^ is a relative correction, i.e. the correction divided by the main amplitude 
(which in our notation is equal to (— a )). The integrand in Eq.(|8|) includes contributions 
from all diagrams, each fi(z), z = k 2 originates from the (i)-th diagram in Fig.l. The above 
described procedure presents the functions fi(z), fiiz) via the two-dimensional integrals 



fi(z) = - du dv[--—z + - + 2ln — ) 
z Jo Jo \D V z u u z / 



h{z) 
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3(1 -w) 



2D 2 D 3 u 3 z 2D 
where the symbols A,B,...F in the integrands are defined as follows 
A = (1 - w)(4 - u + 2v) + ( (1 - it) (2 - v) + 2v(l -v)\z 



(9) 
(10) 

(11) 



B 
C 
D 
E 

F 



(1 - u)v ( (1 - u){2 + u - 2v) + v(l - v)z) 
4- (1 -u)(4-u + 2v) , 
u 2 + v (1 — v)z , 

(1 - u) [3(l + (1 - u) 2 ) + v(l - v)z 



[l-u)v 



(1 - u) 2 (l - 3v) + 1 - v + v 2 {\ - v)z 



The integration over u, v in (|9D, (|TTJ|) arises from the integral in Eq.(|6|) (and a similar expres- 
sion for the diagram (2)), u = £1 +£2, v — £,2- Diagram (3) in Fig.l is expressed directly via 
the mass operator calculated by Feynman pTfl , resulting in 



fs(z) 



2z 



2g(z) + zg'(z) 



2(1 -z) 



z + 



4z-4 



1-z 



\nz) +1 . 



(12) 
(13) 



Eqs.(P),(pH|), and (|T^|) are derived using conventional renormalization procedure to deal with 
the logarithmic ultraviolet divergences. In the Feynman gauge each one of the diagrams in 
Fig.l possesses also the infrared singularity oc In A, A — > 0. It is well known that for problems 
of the type considered the infrared divergence should not manifest itself, and, indeed, we 
find that it cancels out for the sum of all diagrams. This cancellation allows us to derive 
convergent integral representations (||),(|ll]), and to express ( |i~2| ) in terms of a well defined, 
finite function g(z) ([Uf). Since the functions A,B,...F in Eq. (|TT|) are all polynomials, a 
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large number of, probably all, integrations in (PD,(|10|) can be carried out in an analytical 
form. However, our main goal here is a numerical value for the factor C PNC . We rely, 
therefore, on numerical calculations whenever it is more convenient. Calculating integrals 
in Eqs.(|9D,(|lOD we find the functions f%(z) and f^z). [] Combining them with f%{z) from 
(|T2|) we find the function f(z) = Z)f=i fi( z ) m the integrand in Eq.(|J). Its asymptotes are 

f () - - f- m ^+| ? z^O, , , 

3 X \ (In z + const) \ , z — > oo . ^ ' 

Finally, calculating with the found function the integral in Eq. (H) we derive the e-line radia- 
tive correction (|?|) 

C N iPne = -l-97^. (15) 

For heavy atoms the found correction is negative and large. This qualitative result agrees 
with Ref. [pl|. Moreover, numerical results for the Cs atom are also in good agreement. 
Eq.(|i5|) predicts that the e-line correction is —0.6%, which is close to —0.7(2)% found in 
0. 



Eq.flTo]) shows that the e-line correction is large due to a large coefficient ~ 2.0 in its right- 
hand side. Naively one could expect this coefficient to be smaller, of the order of ~ 1/tt. It is 
interesting that similar "numerical enhancement" happens for the e-line radiative correction 
for the energy shift that is due to the finite nuclear size (FNS). Analytically this correction 
was examined by Pachucki [^4j] and Eides and Grotch ||25|| . In our notation their result can 
be written as 

5™f ne = C FNS a 2 Z = 2.978 a 2 Z , (16) 

where C FNS = (3/2) 1.985 = 2.978, the coefficient 1.985 is taken from Eq.(9) of JZ3J. Eq.(|T§ 
shows, that, indeed, the e-line corrections to the FNS energy shift are governed by the 
large coefficient ~ 3.0, similar to Eq. fll~5|) for the corrections to the PNC amplitude. Fig. 2 
examines this similarity in more detail. It shows data available for relative e-line corrections 
for the two problems mentioned above, namely for the PNC amplitude and FNS energy 
shift. The linear in Z approximations (p)|) , (|I6|) (that are valid for sufficiently small values 
of Z) are compared in this figure with results of numerical calculations available for large Z. 
The e-line FNS corrections were calculated in Refs. |f2Tj|-|2"3]] . Specifically, the results for FNS 



shown in Fig.2 were extracted in |18| from |22| (using also Ref. |26j for Z = 55), see details 
in The data shown for the PNC amplitude are taken from the same paper \TE\ that 
explains that the correction to the PNC amplitude is equal to the average of the corrections 
to the FNS energy shifts for Sy2 an d P1/2 levels. The latter ones (i.e. the corrections to 
the FNS energy shifts of Pi/ 2 levels) were also calculated in Refs. p0|-p3|. The results for 



the PNC amplitude shown in Fig.2 are based on data of p2[ (and f26fl), see [18| for details. 



1 A simple way provides an expansion of 1/D in Eqs. (|9|) , (|To|) in powers of [zv 2 /(u 2 + vz)\. For 
each term in this series Mathematica [19| reliably performs integrations in (p|) , (|l0|) in an analytical 
form. The procedure rapidly converges. 



5 



We observe very close quantitative similarity between corrections to PNC and FNS. In 
both cases the linear approximations (|15|) and ([TJ]) predict large negative corrections, which 
agrees with results based on numerical calculations for heavy atoms. Numerical validity of 
the linear approximations seem to be limited by the region below Z = 55, for higher Z they 
underestimate the effect. The Cs atom lies on the border, where results of small-Z and 
large- Z approaches agree reasonably well. 

Numerical data used in Ref. |TE[ incorporates an error that increases for smaller values 



of Z, see Tables III and IV of [ 22]. In order to reduce an impact of this error we can combine 



Eq.(|T5D of the present paper with results of \TE\. Let us approximate the nonlinear terms 
omitted in Eq.(|l~5|) by a simplest quadratic term and choose a corresponding coefficient to 
reproduce the results based on numerical data for very large Z, Z ~ 90, where the numerical 
errors are small. This gives us the following interpolating formula for the corrections to the 
PNC amplitude 

CTne, int = ~l-97a 2 Z(l + 1.55 aZ) . (17) 

Fig. 2 shows that any reasonable interpolation between data available for large-Z and small-Z 
regions would produce a similar pattern. We conclude therefore that Eq.flTTD should give 
the most reliable numerical data. For the Cs atom we find from it the correction —0.9(2)%, 
where we adopt the error 0.2% of |1| (compare this result with —0.6% of Eq.([L5|) and 



—0.7(2)% of Ref. |18] mentioned above) 



The standard model value for the nuclear weak charge for Cs |28[ is 

Q w ( 133 Cs) = -73.09 ± (0.03) . (18) 

Ref. H refined previous calculations of Ref. || extracting from the experimental PNC am- 
plitude of Ref. || the weak charge 

Q w ( 133 Cs) = -72.18 ± (0.29) cxpt ± (0.36) theor , (19) 

with the theoretical error 0.5%. It is consistent with Q w ( 133 Cs) = —72.21 ± (0.28) expt ± 
(0.34) theor that was adopted in |Tj| by taking the average of the results of Refs. [f§-f|, and 
accepting the theoretical error 0.4% of 0. The weak charge in Eq. (|i~9"D deviates from the 
standard model ( |I~8"D by 2.0a. Taking from Eq.([l7]) the value —0.9(2)% for the e-line radiative 
correction, we derive from Eq.(|l9"D 

Q Ty ( 133 Cs) = -72.83 ± (0.29) cxpt ± (0.39) theor , (20) 



which brings the experimental results of within the limits of the standard model (|T8|). 



(Note that even the smaller value for this correction —0.6% that follows from the linear 
approximation Eq. (|T5"D justifies the latter conclusion.) 

The calculations presented show that the QED self-energy and vertex radiative correc- 
tions bring the experimental results of Wood et al [^] within the limits of the standard 
model. 

This work was supported by the Australian Research Council. I wish to thank Victor 
Flambaum for numerous discussions, and K.T.Cheng for Ref. 
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APPENDIX A: 



The calculations performed use the following algebraic relations 



(75) = , (7^75) = , (7^75) = -(p»a» ~ Vvdp) 
(7m7^7a75) = -{a^gux - a u g^ x + a x g^) - it^ u \ p a p , 

where the (X) = Tr(Xp), and conventional integrals 



(Al) 



/ 



(q 2 + 2qs - M 2 ) n 
d 4 q 

(q 2 + 2qs - M 2 ) n 

q^q" d A q 
(q 2 + 2qs - M 2 ) 3 
qVq"q x d 4 q 



{-l) n m 2 c n 
(M 2 + s 2 ) n - 2 



-l) n - l m 2 c n s» 
(M 2 + s 2 ) n - 2 ' 



n > 2 , 
n > 2 



-in 



n u 



s^s 



1 



M 2 + s 2 
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In 



(-l)^7r 2 c„ 

(q 2 + 2qs - M 2 ) n (M 2 + s 2 ) 



n— 3 



s fl s' / s x 
'M 2 + s 2 



A 2 _3 N 

~ 2 + s 2 ~ 2 j 

fiu „A 1 n fJ,\ Q v 



+ 



2(n-3) 



(A2) 



n > 3, 



where c„ = l/[(n — 1)(^ — 2)], M is a scalar, s is a four-vector, gs = (qs), and A is the 
ultraviolet cut-off. 



7 



REFERENCES 



M.A.Bouchiat and C.Bouchiat J.Phys. (Paris) 35, 899 (1974); 36, 493 (1974). 
S.L.Gilbert and C.E.Wieman, Phys.Rev. A 34, 792 (1986). 

C.S.Wood, S.C.Bennett, D.Cho, B.P.Masterson, J.L.Roberts, C.E. Tanner, and 
C.E.Wieman, Science 275, 1759 (1997). 

V.A.Dzuba, V.V.Flambaum, and O.P.Sushkov, Phys. Lett A 141, 147 (1989). 
S.A.Blundell, J.Sapirstein, and W.R.Johnson, Phys. Rev. D 45, 1602 (1992). 
M.G.Kozlov, S.G.Porsev, and I.I.Tupitsyn, Phys. Rev. Lett. 86, 3260 (2001). 
V.A.Dzuba, V.V.Flambaum, and J.S.M. Ginges, |hep-ph/01110T9l (2001). 



V.A.Dzuba, V.V.Flambaum, and J.S.M. Ginges, |hep-ph/0204T34] (2002) 



S.C.Bennett and C.E.Wieman, Phys. Rev. Lett. 82, 2484 (1999); 82, 4153 (1999); 83, 
889 (1999). 

A.Derevianko, Phys. Rev. Lett. 85, 1618 (2000). 

V.A.Dzuba, C.Harabati, W.R.Johnson, and M.S.Safronova, Phys. Rev A 63, 044103 
(2001). 

O.P.Sushkov, Phys. Rev. A, 63, 042504 (2001). 

W.R.Johnson, I.Bednyakov, and G.Soff, Phys. Rev. Lett 87, 233001-1 (2001). 
A.I.Milstein and O.P.Sushkov, |hep-ph/0109257| (2001). 



M.Yu.Kuchiev and V.V.Flambaum, |hep-ph/02050T2] (2002). 
W.J.Marciano and A.Sirlin, Phys. Rev. D 27, 552 (1983). 
B.W.Lynn and P.G.H.Sandars, J. Phys. B 27, 1469 (1994). 
M.Yu.Kuchiev and V.V.Flambaum, [hep-ph/ 02061241 (2002). 

S. Wolfram, Mathematica. A system for doing mathematics by computer. (Reading, Mass. 
Addison- Wesley Pub. 1991). 

W.R.Johnson and G.Soff, At. Data Nuc. Data Tables 33, 405 (1985). 
S.A.Blundell, Phys. Rev. A 46, 3762 (1992). 

K.T.Cheng, W.R.Johnson and J.Sapirstein, Phys. Rev A 47, 1817 (1993). 

I.Lindgren, H.Persson, S.Salomonson, and A.Ynnerman, Phys. Rev. A 47, 4555 (1993). 

K.Pachucki, Phys. Rev. A 48, 120 (1993). 

M.I.Eides and H.Grotch, Phys. Rev. A 56, R2507 (1997). 

K.T.Cheng. Private communication (2002). 

R.P.Feynman, Phys.Rev. 76, 769 (1949). 

D.E. Groom et al, Eur. Phys. J. C 15, 1 (2000). 



8 



FIGURES 



FIG. 1. The QED self-energy and vertex corrections, called e-line corrections in the text, to 
the PNC matrix element. For each diagram one of the wavy legs shows the Coulomb interaction 
with the nucleus, another one the weak PNC interaction with the nucleus. Each diagram represents 
all possible Feynman diagramms with the given topological structure. 



FIG. 2. The thin dotted line, thick dotted line, and solid line: e-line corrections to the 
PNC amplitude that follow from Eq.(15), Ref.[18], and Eq.(|T7|). The thin dashed line and thick 
dashed line: the e-line corrections to the FNS energy shifts that follow from Eq.(16) (derived 
from [25]), and Ref.[18] (based on calculations of [22]), the dashed-dotted line: the interpolation 
C-iL int = -2.978 a 2 Z(l + 0.85 aZ) between the two latter lines. 
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Fig.2 



